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^SJ . Calculating the frequency dependent dielectric function for strongly coupled plasmas, the rela- 

■ tions within kinetic theory and linear response theory are derived and discussed in comparison. In 

' this context, we give a proof that the Kohler variational principle can be extended to arbitrary 



I. INTRODUCTION 



frequencies. It is shown to be a special case of the Zubarev method for the construction of a non- 
equilibrium statistical operator from the principle of the extremum of entropy production. Within 
\l ■ kinetic theory, the commonly used energy dependent relaxation time approach is strictly valid only 

<n: for the Lorentz plasma in the static case. It is compared with the result from linear response 

theory that includes electron-electron interactions and applies for arbitrary frequencies, including 
(— I ■ bremsstrahlung emission. It is shown how a general approach to linear response encompasses the 

^ I different approximations and opens options for systematic improvements. 

PACS numbers: 52.25.Dg,52.25.Fi,52.25.Mq,52.27.Gr 

Q ■ Different approaches have been elaborated to evaluate the response of a plasma to external time and space dependent 
, electric fields. This applies, e.g., to absorption and emission of radiation (in particular bremsstrahlung), Thomson 
scattering, and stopping power. The dielectric function e{k,uj) depending on the wave number k and frequency uj as 
i-S^ . the central quantity is related to the polarization function, the dynamical conductivity or the dynamical structure 
factor. The random phase approximation (RPA) is improved if collisions are taken into account. In this context, a 
non-local dynamical collision frequency is introduced [H-Q • Alternatively, the concept of a local- field factor Q can be 
CO ■ extended to dynamical local- field corrections 0, [^-Q . In the present work, we focus on the dynamical conductivity and 
^ ' restrict ourselves to the long-wavelength limit fc — >■ 0, i.e., the response of a charged particle system to a homogeneous, 
time-dependent electrical field. 

A well-known expression for the dc conductivity of a fully ionized plasma in the classical, low- density limit has been 
given by Spitzer and Harm [9] within kinetic theory (KT). Further approaches by Lee and More Stygar jTlj . 
and others improved the electron-ion interaction using the relaxation time approach. However, to recover the Spitzer 
result for the conductivity, electron-electron collisions have to be taken into account. This is not consistently possible 
within the relaxation time approach [l^] . We discuss a general approach that allows also for a systematic treatment 
of electron-electron collisions. 

The investigation of time-dependent fields is somehow difficult in KT, too. Often, a combination of the collisionless 
_ kinetic equation with the relaxation time ansatz is used, see Landau and Lifshits Dharma-wardana or 
. Kurilenko et al. [l|, [l^ . It has been emphasized by Landau and Lifshits f\3\ that such an approach is only applicable 
\^ in the low frequency limit. The high-frequency region, where bremsstrahlung is relevant, has to be treated in another 
way. In this work, we present general expressions applicable to arbitrary frequencies of the external field. 

In linear response theory (LRT) , the Kubo formula [l^ was considered as a promising approach to the dynamical 
conductivity in dense, strongly interacting systems at arbitrary degeneracy. A generalized approach to non-equilibrium 
processes has then been given by Zubarev et al. [iTj that will be applied here. It relates transport properties to 
equilibrium correlation functions such as current-current or force-force correlation functions. Different methods can 
be applied to evaluate these correlation functions such as numerical simulations, density functional approaches [13, [ill > 
or analytical expressions derived from perturbation theory |l9l - l2ll| . Note that also strict results such as sum rules can 
be employed to construct the dynamical structure factor, see [22l424l |. We will show how consistent approximations 
are obtained from a general scheme of non-equilibrium statistical physics and systematic improvements can be given. 

In the present work, we will restrict ourselves to homogeneous systems and therefore do not consider any dependence 
on the position r in space, e.g. due to external potentials, in addition to the homogeneous, time-dependent electrical 
field that is treated as perturbation. The focus is on the generalization of relations which were originally derived in KT, 
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see Sec. HIl Starting from LRT, see Sec. IIIIl a generalized Boltzmann equation with a frequency dependent collision 
term is derived. In Sec. llVi a variational approach is applied for the solution of the generalized linear Boltzmann 
equation. Similar to the use of polynomials [ll, [1^ [l^ to solve the static Boltzmann equation, we consider moments 
of the single-particle distribution function to find approximate solutions. Furthermore in Sec. Sec. |lTl different 
limiting cases such as the dc conductivity and the high-frequency limit of the absorption coefficient are considered. 
The dynamical conductivity from KT using an energy dependent relaxation time that has often been used in the 
literature, is compared with the rigorous treatment within LRT. Conclusions are drawn in Sec. [V] 



II. KINETIC EQUATIONS 

A. Single-particle distribution function 

We consider neutral Coulomb systems that consist of charged particles such as electrons and ions. Response to an 
electromagnetic field is described by the dielectric function, taken in the long-wavelength limit here, 

lim e(fc,w) = 1 + — cr(w) , (1) 

fc— s-oo £qUJ 

or the dynamical conductivity (t(w). Treating the Coulomb interaction in mean-field approximation, the random 
phase approximation (RPA) is obtained. To improve RPA, one has to include collisions. A standard way to treat 
collisions is the Boltzmann equation where the interaction between the constituents leads to the collision term. As a 
consequence, a dynamical collision frequency can be introduced according to a generalized Drude formula 

= (2) 



with the plasmon frequency — ^ e^nj (com), where n is the electron density and m the (reduced) mass. The 
collision frequency v(lo) should be a complex, frequency dependent quantity in order to satisfy sum rules of the 
dieletctric function. This is in contrast to a static relaxation time r — l/i/, as used in the kinetic approach and will 
be explained in more detail below. 

Taking the adiabatic approximation, N electrons interact with singly charged heavy ions that are considered as 
external potential. The Hamiltonian with the electronic degrees of freedom only, is 

= ^ Ephlap + ^ V^^{q)a\^^a.p + ^ ^co('?)4i+9°P2-«"P2api (3) 
V pq pipiq 

with Ep — fi'^ / {2m) . The interaction with the ions 14i((j) = — ^(g) 6xp[i(f • Rj] describes Coulomb potentials 
V{q) — I (eoilog^) due to various ion sites Rj what leads to the structure factor S{q) — j exp[i(f-(^i— ^j)]). 

VLq is the normalization volume. The electron-electron interaction is given by the Coulomb interaction Vcdq) — y{q)- 
The account of the ion dynamics is straightforwardly taken into account within a two-component plasma [3| , but the 
notations become more complex and will not be given here. 

For the derivation of kinetic equations, in particular the Boltzmann equation, we consider the electron single- 
particle distribution function /(p, i) = Tr{np/5(i)} = {fipY that is the quantum statistical average, taken with the 
non-equilibrium statistical operator of the single-particle occupation number operator hp = a^dp of momentum 
Tip. Considering homogeneous systems, the density matrix is diagonal with respect to the wave vector p. Spin variables 
are not explicitly given unless it is pointed out. Subsequently, the single-particle distribution function does not depend 
on the position r either. 

In thermal equilibrium, the single-particle distribution function /o(p) — Tr{np/5o} is calculated with the grand 
canconical statistical operator po = exp[— — p/V)]/Tr{exp[— /3(i? — ^JiN)]}■ Neglecting the interaction term, we 
find the ideal Fermi gas with distribution fp — [exp(/3(£'p — ^)) + 1]^^. Under the infiuence of an external perturbation 
iJ*xt: the single-particle distribution function f{p,t) is modified. Its deviation 

Sfip, t) = /(p, t) - Up) = Tr {5hp p{t)} (4) 

from the equilibrium distribution /o(p) is the average of the fiuctuations of the single-particle occupation number 
Sfip = hp — fo{p). The time dependence of the single-particle distribution function /{p^t) is determined by the 
non-equilibrium statististical operator p{t) as shown in the following section. 
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Alternatively, the dynamics of the single-particle distribution function can be determined from a hierarchy of 
equations of motions for the many-particle distribution functions. Truncating the hierarchy, a kinetic equation [iSj is 
obtained with the following structure 

§-^f{p,t)^B[fip,t)]+C[fip,t)] (5) 

describing drift in the single-particle phase space via drift term D[/(p, t)], and collisions that are caused by the 
interaction between the particles. The collision term C [f{p,t)] is related to higher order distribution functions due to 
the interaction mechanisms within the system. To obtain closed kinetic equations, the higher distribution functions 
are expressed in terms of /(p, t). 

In the following, we consider a homogeneous system under the influence of an external time dependent electric field 
E{t). The total Hamiltonian H^^^ = H + i?*xt contains the interaction with the external field H^^^ = —eE{t) ■ 
for the electron position operators r^. From the respective external force eE(t), the drift term follows as 



D [fiP: t)] - ■ ^_f{p, t) « - /,) E{t) ■ p (6) 

n op m 



in first order with respect to the external field E{t), with (3 — l/{k^T). Expressions for the collision term C [f{p,t)] 
will be given below. 

With the distribution function /(p, t), the current density is given by 

Jit) ^^Y.^pf{p,t) = ^Pi(i). (7) 

p 

The total momentum Pi{t) is the first moment of the distribution function. In the following, we also consider the 
operators of arbitrary moments 

P. ^Y.^PEmpf-^y^fip (8) 
p 

where pe = P ■ E /\E\ denotes the component of p in the direction of E. 

Arbitrary time dependence of an electric field can be expressed by superposition of harmonic time dependences. 

Within linear response, each component E{t) = ■^E{uj) exp(— iwt) -I- c.c. causes an induced single-particle distribution 
function 

^fiP, t) = ^Sf{p,uj) exp{-iujt) + c.c. (9) 
with the same time dependence. The dynamical conductivity follows from j{io) — a{uj)E as 

a{u:) = ^^Y.^P^ W'^) ■ (10) 
mE "0 — 

p 

Note that all Fourier components marked with tilde, e.g. Fp^ are frequency dependent in general. The dependence 
on uj will be omitted in some of the following expressions for them to be more compact. 



B. Relaxation time approximation and dynamical conductivity 

To start with an analytically solvable example, we first discuss the solution of the kinetic equation ([S]) for the 
Lorentz model where the electron-electron interaction in the Hamiltonian ([3]) is neglected. Considering a constant 
electric field, the distribution function f{p,t) = f{p) is static. In the standard treatment, see [27|], the collision term 
reads 

CLorentz [f{p}] = ^ {f{p)w,i{p,p') [1 - f (p)] - f {p)w,,{p' ,p) [1 - /(p")]} • (H) 
pi 

The transition rates can be determined in Born approximation from the golden rule, Wei{p,p') = {2Tr/h)\Vci{\p — 
p'\)\'^S{Ep — Epi). Since the energy of electrons is conserved in adiabatic approximation, a relaxation time Tp is 
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introduced via an ansatz for the linear term of the expansion of the distribution function f{'p) = fp — Fpj-^rfp- 
analogy to the drift term ([6]) we assume 

Sfip) - -^TpE-p, (12) 
m 

which realizes the linearity with respect to the external field E. For isotropic systems, Tp is a scalar depending only 
on the modulus of p. Inserting Eq. ()12|) into the collision term ([Tl]) and taking into account the detailed balance 
in equilibrium Wci{p,p')fp'{l- — fp) = Wci{p',p)fp{l — fp') as well as the energy balance of the transiton rates, the 
collision term (|lip is 

CLorcntz [fip]] - - ^ W,-Mp')fpil " fp){Fp - Fp,) ^ -5f{p)/Tp (13) 
p' 

For the kinetic equation ([5]) with the drift term Eq. ([6]) we then find 

E-p = - ^Wci{p,p')-^ E ■ {Tp>p' -Tpp) = -Tp^Wc\{p,p + q)E ■ q (14) 

with q = p' — p. With the golden rule for the transition rates given above and S{q) « 1, |t4i(<z)P ~ NV'^{q), the 
energy dependent relaxation time can be calculated 

Tp n ^ E-p 

The q integral in Eq. (jlSp can be performed using spherical coordinates where p is in z direction, E in the x — z plane. 
It is convergent only in the case of a screened Coulomb potential. Using the statically screened Debye potential 

e2 

''^0(9) = o ^ 2 I — 2T' '*D = /3"eVeo, (16) 

we find the energy dependent collision frequency 

-1 e'' m / n r I b \ 

Vr, = T„ = n — in V 1 + 17 

with b = 4p^/k^ in the Coulomb logarithm. The static conductivity is determined from Eq. ([TOl) . cj = 0, as 

e^n^ . 1 



*-^dc,Lorcntz 



m 



2 P-FT X! Pe fp{^ - fp) = eoWplTLorontz = • (18) 

^ ^0 ^ ^^Lorentz 



P 



We introduce the average relaxation time TLorcntz and the static collision frequency t'Lorontz = l/^Xorcntz- 

We are now interested in extending the static case, Eq. ([18]), by evaluating the permittivity e(aj), Eq. ([T]), or the 
dynamical conductivity, Eq. (jlOp . From the kinetic equation ([5]) with the drift term Eq. ([6]) we derive the frequency 
dependent Boltzmann equation 

eh ^ r ~ 1 

-iujSf{p,Uj) = /^^^(w) •p/p(l - /p) + CLorcntz S f {p, Oj) . (19) 

TO i- f |_ J 

In a standard approach, see, e.g.. Landau and Lifshits [27| . it is proposed to extend the static case to the dynamic 
case assuming that the relaxation time is the same as in the static case, see Eq. (1131) . Subsequently, the following 
relation is derived, 

- -) SRp, io) - —l3E{u;) ■ pfp{l - fp) , (20) 

TpJ TO 

SO that for the dynamical conductivity ([T0|) follows (spin factor 2, p^ — > p^/3 for isotropic systems) 

2e'n'l3 f d'p p'fp{l-fp) 
^kt(^) = ^—^j (2^— • ^^^^ 
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This result can be interpreted as a Vlassov approach where the frequency uj is replaced by a complex frequency 
w + i/Tp. However, the introduction of an energy dependent, static relaxation time is an approximation that cannot 
be applied, in particular, at high frequencies, where bremsstrahlung emission is expected. Note that it is not possible 
to give an explicit expression for a frequency dependent collision frequency as desired for a generalized Drude formula 
according to Eq. ([2]). Furthermore, inelastic collisions such as electron-electron interactions are not taken into account 
by a collision time ansatz. Further evaluation of Eq. ()2ip is given in Appendix |E1 results are shown in Fig. [2] and 
discussed below. 



III. LINEAR RESPONSE EQUATIONS 
A. Linear response theory 

To evaluate the response Q to an external perturbation H^-^^ , we determine the non-equilibrium statistical operator 
p{t) within a generalized linear response theory. The conceptional ideas and main expressions relevant for the further 
analysis of the single-particle distribution function will be given here according to 0, [28l - [30| . 

We introduce the relevant statistical operator 

= ^^e-^(^'^*)+^"^"W^", =Tr|e-^(^-^^)+^-.^"(*)^4 , (22) 

as a generalized Gibbs ensemble which is derived from the principle of maximum of the entropy 

S{t) = -A:BTr{p,ei(t) ln[p,ei(t)]} , (23) 
where the Lagrange parameters /?, /i, F^it), which are real valued numbers, are introduced to fix the given averages 

Tr {B^m} = =Tr{B„p,ei(i)} • (24) 

These self-consistency conditions mean that the observed averages {BnY are correctly reproduced by the hermitean 
Prciit)- Similar relations are used in equilibrium to eliminate the Lagrange parameters /3 and fi. In linear response, 
the response parameters f „ (t) are considered to be small so that we can solve the implicit relation (j24|) expanding up 
to first order. 



Prc\{t) = 



1 + V F„{t) f dAe-'3^(^-^^)<5B„e''^(^-^^) 



Po ■ (25) 



Note that the expansion of .^roi(0 i'^ Eq. leads to the subtraction of the equilibrium average in 5Bn = Bn — (-Bii)o- 
The average fluctuations can now be explicitly calculated by inserting Eq. (|25p in Eq. (IM)) . 

{SB^y = ^((5S„, 6Brn)Frait) , (26) 
rn 

where we introduced the Kubo scalar product 

(i, B)^ dX Tr jist {inf3X)po } . (27) 
The time dependence A{t) = Q^Ht/h^^-iHt/h 

is given by the Heisenberg picture with respect to the system Hamiltonian 

H, and A = i[H,A]/n. 

A statistical operator for the non-equilibrium is constructed with the help of the relevant statistical operator (|22p , 
see App. [X] Expanding up to the first order with respect to the external field E and the response parameters Fn, 
where F„(i) = Re{F„(a;)e~*"*}, we arrive at the response equations 
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with z = uj + ie, the total momentum of electrons P — 'Y^^Tipfip, and the Laplace transform of the correlation 
functions, 

{A]B^ ^ dt e*^* = ^ di e*^*^ dA Tr{i(t-i;i;3A)#po}- (29) 

Considering Nb relevant observables i?„, Eq. (|28|) is a system of Nb linear equations to determine the response 
parameters Fn for a given external field E. It is the most general form of LRT, allowing for arbitrary choice of 
relevant observables i?„ and corresponding response parameters Fn- We show below that, with respect to kinetic 
theory, the first two terms on the left hand side of Eq. ([28| can be identified as a collision term, while the right hand 
side represents the drift term due to the external perturbing field. 



B. Generalized linear Boltzmann equations 



In kinetic theory, the non-equilibrium state is characterized by the single-particle distribution function f{p,t). In 
order to derive expressions in parallel to the kinetic theory, we choose the fluctuations Sfip of the single-particle 
occupation number, see Eq. as relevant observables Bn- The modification of the single-particle distribution 
function can then be calculated straight forwardly according to Eq. (|26|) 

(30) 



Tr {/5rei(i) Sfip} = ^ {Shp, Sfip,) Fp, (t) = df{p, t) . 



The Lagrange multipliers Fp(t) — Fp{uj) exp(— iwt)/2 -I- c.c. are determined from the response equations Eq. ([28l 
arrive at the generalized linear Boltzmann equations {Shp — fip) 

eh 



We 



E 

p' 



\^(5np,np,) + 



^{Sfip,6hp,) + (^np;6hp,^ | Fp> = — /3 ^{Sfip, hp„) + (hp] n. 

p" 



-II 
P 



E. (31) 



The time derivative of the position operator in -ff leads to the total momentum hpi — m and subsequently 
to the right hand side of Eq. (pij) . We analyse the different terms of Eq. ([5T]) below and compare with the kinetic 
equation Eq. ([5|), considering the Born approximation. Notice that this result can be extended by introducing 
stochastic forces [2^ if we go beyond the Born approximation. Further relevant observables beyond the single-particle 
occupation numbers can be included in order to characterize the non-equilibrium state, such as long-living correlations 
and formation of bound states. It is possible to go beyond the Boltzmann equation if higher correlations such as bound 
state formation are included into the set of relevant observables. 
We give the entropy as obtained from Eq. ([23)) 



S{t) 



- ln[Zrci(t)] - f3{H -^iN)+Y, Fp{t) 



5o(/3, A^) - fcs ^ Fp{t) Sfip, t) (32) 



in first order oi Fp{t). The entropy in the thermodynamic equilibrium is denoted by So{l3,fi). With Eq. (1301) we find 
that the entropy decreases in non-equilibrium because SS{t) = — X]pp' ^p' (0 (Snpi , drip) Fp{t) < 0. The proof is given 

using the spectral density for F{t) = Fp{t)Shp, see [iTj]. With the eigenstates {H — iiN)\n) — En\n) of the system 
Hamiltonian we have 



5S{t) 



F{t),F{t) 



1 ^e-I^E- 



Zol3 



'i3E„ 



En — E„i 



-\{n\F{t)\m)\'<0. 



(33) 



This result corresponds to the second law of thermodynamics that the entropy of the many-particle system exhibits 
its maximum in the equilibrium state. 



C. Evaluation of equilibrium correlation functions, Born approximation 



Quantum statistics provide us with different methods to calculate correlation functions in thermal equilibrium such 
as perturbation theory and diagram techniques. Applying perturbation theory with respect to the interaction. Wick's 
theorem can be used. We find in lowest order for the Kubo scalar product, Eq. (P7)) . 

{■hp, hp') = Tr |po ftp' V«p'^p} = /p'/p + /p(l ^ /p)'^pp' (34) 



7 



so that {5rip,5hpi) = {5hp,hpi) = fp{l — fp)5ppi. The remaining Kubo scalar product vanishes, {5hp,hp') = 0, as 
shown from the Kubo identity Eq. (jA3| with C = Sfip, and ([np',ftp])o = after cychc invariance of the trace. 

For the deviation of the single-particle occupation numbers from equilibrium we find from Eq. pOI) that df{p,t) = 
Fp{t)fp{\ — fp) which is equivalent to the expansion (fT2|) in kinetic theory. Thus, we solved the self-consistency 
condition ([24]) to eliminate the Lagrange parameters Fpit). According to ([9|), the Fourier components 

5}[p,u:) = fp[l- fp)Fp{u) (35) 

are complex amplitudes, containing in general a phase factor. 

The equation of motion that leads to the generalized linear Boltzmann equation pip allows to relate the response 
to the external field. The right-hand side is the drift term that contains the external field. In Born approximation, we 

can neglect the correlation function (^np\ fip"'^ because it is of higher order of interaction compared with (Sfip, fip"). 

Then, the right-hand side of Eq. (|31l) reads 

Bp^^Pfp{l-fp)p-E (36) 

in agreement with Eq. By the same argument we have the term due to the explicit time dependence 

— |^(5np, (5np') + l^fip^Shpi^ Fp' = —iujSf{p,u!) = —iVtpFp , (37) 

p' 

with Vlp — ^fp{^ — fp)- Note that the correlation function (hp;hp") is eliminated introducing stochastic forces 



[ITl |29j so that the result —iujdflp, uj) holds also beyond the Born approximation. 
The remaining term in Eq. pip describes the collision integral. 



Cp ^ -'^(hp;hp>\ Fp> = -'^Cpp'{uj)Fp'. (38) 



P' 



It is evaluated in Born approximation, see Appendix iBl with the generalized Onsager coefficients Cpp>{uj) — (w) 
Cpp,{u}), leading to 



+i nr+i 



{Ep — Ep+q)/fl UJ — {Ep — Ep+q)/fl 



tt5[lo + ^{Ep - Ep+q)^ +TrS ^ - ^'■"^^ ~ ^p+'i 

[Sp' ,p+q ~ Sp' ,p] J (39) 



'fpfpii^ fpi-q fp+q) fp+qfpi-q{^ fpi fp) 



h /3{Ep+q + Ep^^q - Ep^ - Ep) 

{^p' .p+q ^P\Pl ^P' rPl ^p' -p] 5 (^^) 



Pi, 9 

i i 



UJ + ie + Ap^p^^q UJ + ie - Ap^p^^q_ 

where Ap^p^^q = {Ep^q + Ep^_q — Ep^ — Ep)/h. Exchange contributions have been discarded, see Appendix [B] The 
decomposition of £pp, (w) in real and imaginary part is analoguous to C^p,{uj). 

In conclusion, the generalized linerized Boltzmann equation pip can be given in the same way as assumed in the 
relaxation time approach, see Eq. (|19p . 



- lUjSfip, uj) = ^/5/p(l -fp)p-E~J2 ^PP' (^)^P' = Dp + Cp[Sfip, uj)] (41) 

p' 

with Eq. ((37|) and the drift term (l36t . after replacing the response parameters Fp in the collision term (|38| by the 
single-particle distribution according to Eq. psp . This holds for arbitrary frequencies uj and degeneracy, see Appendix 
A. At zero frequency, the collision integral ([TT|) of the Lorentz plasma is recovered if calculations are taken in Born 
approximation and restricted to the electron-ion interaction only. At arbitrary frequencies, the collision integral 
becomes a complex quantity in contrast to the scalar relaxation time. Real and imaginary part are connected via 
Kramers-Kronig relations. The Born approximation can be improved in a systematic way if the correlation functions 
are evaluated in higher orders with respect to the interaction. A Kubo-Greenwood formula can be derived that 
expresses the collision term by T matrices [1^, H^l • 



8 



IV. SOLUTION OF THE GENERALIZED LINEAR BOLTZMANN EQUATION 

A. Variational principle 



Having derived explicit expression for the Onsager coefficients Cpp' in Born approximation, Eqs. p9p . (PU]) . we can 
now determine the response parameters by solving the generalized linear Boltzmann equation (j4ip given as 



iilpFp{uj) = Dp - £pp> {u)Fp> {u 



(42) 



As a further constraint on the response parameters Fp, we consider the entropy leading to a variational problem as 
follows. 

We determine the time derivative of the entropy, Eq. (l32t . The time dependent term reads 



^p[E] (e— * +e-*) -E£pp,(a.)f'p,e— * -5]/:;,,(a.)i?;,e^' 



2iujt ^—2iujt 



(43) 



arise that 



if we insert the Boltzmann equation (|4T|) for —iuj6f{p) for the last line. Oscillating terms cx e' 
disappear in the time average. The remaining terms cancel, which can be directly seen, if replacing Sf{p) by the 
Langrange multipliers Fp using Eq. (|35|) . Thus the total entropy is constant in the average over a period of time, 
dS(t)/dt — 0. However, even in the time average, there is an entropy production which is dissipated as entropy export 
due to the external field in the drift term. We have 



Therefore, let us consider the functional 



pp' 



Sint[Gp] = ^^GpCpp>{uj)Gp' + c.c. 
pp' 



(44) 



(45) 



for any function Gp that obeys the constraint 



{^)Gp. 



= 



(46) 



that can be considered as an integral over the Boltzmann equation (|42p . It is easily shown that the time averaged 
change of entropy Eq. (j43p vanishes for arbitrary functions Gp that obey the constraint (|46p . The maximum of the 
functional S'intiG'p] occurs at Gp = Fp which is the solution of the linear Boltzmann equation (|42l) . see App. [C]for the 
proof. 

This is a generalization of the Kohler variational principle fsil. l32j| for arbitrary frequencies w. It can be related to 
the principle of extremum of entropy production given by Prigogine and Glansdorff 33]. The static case w = has 
been considered in Refs. [13, Isil. f32l. l34| . Some attempts to extent this to arbitrary frequencies can be found in [35| , 
but, to our knowledge, a consistent approach has not been given until now. 

In order to apply the variational principle given here, one can consider a class of trial functions G'^^''\^v\'p) — 
'Yliu'Lx ^v9v^V) with respect to an arbitrary but finite (iV^) set of linear independent functions guij>)- Determining the 
extremum of iSmtl^i/] leads to an optimal set of parameters — F^^"^ . The extension of the class of trial functions 
to an infinite number of functions then gives the exact result Fp = limAr^_j.oo ^Yl!l=\ ^v^"^ 9v(,v)- 

Alternatively, the relevant observables hp are replaced by a reduced set of A^^ relevant observables — X)p gv{v)hp. 
The solution of the finite system of linear equations (|A4[) then gives the Lagrange multipliers Fy , that can be expressed 
in terms of determinants. This leads to identical results as for the variational principle. In previous papers we used 
a finite number of moments g^ip) = JiPEWEp)^'^~^^^^ according to the general moments ([5]). An alternative basis set 
would be the Sonine polvnomia ls 12511 t hat are appropriate in the static, nondegenerate limit. It has been shown that, 
within perturbation expansion [36|, |37[, results are converging with an increasing number of moments used. 
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B. One-moment Born approximation 



In lowest approximation, we choose with Gp = Figi{p) — FihpE the first moment of the distribution function (|8]) 
as trial function. The variational parameter Fi is fixed by the auxiliary condition (j46p where we insert Eq. (I36p and 
flp from Eq. (j37p . and we find 



eh 



F^hpE—Pfpil - fp)PEE = -luY^iF^fipE? fp{l ~fp)~J2 Fi^PEC%,{oj)F,hp'E . 



(47) 



The electron-electron collisions do not contribute in the one-moment approach because of conservation of total mo- 
mentum. We assume the general structure of the variational parameter 



e/3 



1 



E. 



(48) 



After some calculations given in Appendix [Dl we find the collision frequency for the case of the statically screened 
Coulomb potential Eq. (IT6)) . and S{q) « 1, 



1^0(1^) = «.9do 



dy 



with 



(y^ + n)' 



1 em 



dx 



X w + le — X 



■In 



1 

1 + e^(^/2'+!^)^+'3*' 



48^4 ^27^3 ' 



w 



8m 



(49) 



(50) 



which is valid for any degeneracy. In the non-degenerate limit <C 1, we can expand the logarithm. With e^^ — 
n(27r/3?i^/TO)'^/^/2 — nk^ /2 and spin factor 2, we find 



^'d(w) ^ign dy 



V 



(j/2 + nf- 



dx 



1 - e-'^'^y 

xy{w — xy + ie) 



(51) 



with g = 5degen/2. 

The dynamical conductivity (jlOp can now be calculated with Eq. p5p and the optimized Lagrange parameter 
Eq. (05]) so that Fp ^ FiUpE- We find 



iF "0 — 



(52) 



For isotropic systems, the sum is evaluated as J^pi^Ps)^ fpi^ ^ fp) = Nm/l3, see App. [D] Inserting the derived 
expression (|48| we obtain a generalized Drude type expression, Eq. ([2]), 



o-d(w) 



—iuj + fD(i^) 



(53) 



for the dynamical conductivity. The comparison with (Tkt, Eq. (1211) . will be performed in the following Section. 

It is instructive to investigate the alternative approach where only moments of the distribution function P^, Eq. ([8|), 
are taken as relevant observables _B„, instead of the fiuctuations Sfip of the single-particle occupation operator as 
originally introduced in Subsec. IIIIBI Taking the component of the total momentum of the electrons Pi = hpETip 
in the direction of E as an one-moment approach, we have with Eq. ([T]), (I26p 



{Pi)Fi 



■{Pi.Pi)Fi 



(54) 



The generalized linear Boltzmann equation (|3ip is now reduced to a single equation that reads in Born approximation 
{{Pi)o = in thermal equilibrium) 



Pi; Pi) -zw(Pi,Pi) 

LJ-\-ie 



Pi = {Pi,P,)-i3E 
m 



(55) 
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containing force- force correlation functions as the coUison term. With (Pi, Pi) = Nm/fS, see App. |D] and the 
staticahy screened interaction Eq. (|16p . the expressions for the dynamical conductivity, Eq. (j53p . and the corresponding 
dynamical collision frequency 

^lf''M = ^7r(A;Aw (56) 

is obtained that coincides with the results Eq. (|49p and Eq. ([STj) given above. This is a preliminary result of the 
LRT based on the one- moment Born approximation. Going beyond the Born approximation, we denote {/'^^^^(a;) — 

/3/(m iV)(Pi; Pi)ij+ie as collision frequency of the one-moment approach. Systematic treatments of the perturbation 
expansions are performed with the help of Green's function techniques. In particular, the Gould-DeWitt approximation 
for v^^^^uj) has been performed that accounts for the correction of long-range interaction by dynamical screening and 
considers strong collisions at short ranges [1, [s^] ■ 



C. Higher moment approaches 



An improvement of the dynamical conductivity (j53p can be achieved by extending the set of trial functions or relevant 
observables within the variational approach or the relevant statistical operator, respectively. Using higher order 
moments P^, Eq. ([5]), of the distribution function, converging expressions are obtained for the transport coefficients 
[13, HI]. In particular, higher moments are needed in order to take into account electron-electron collisions. Taking 
higher order moments into account, the change of the dynamical conductivity can be represented by a complex 
function r{u}) so that iy{uj) = r(cj)z^(^i)(a;) 0,11 Hi], 



(57) 



As a special case, we discuss the two-moment approach with Pi, P3 as relevant observables (i.e. particle current 
and energy current). The account of these two functions in p space allows for a better variational approach to the 
single-particle distribution function. For the electrical current density we have with Eq. ([71), (|26p 



(A) - |(A,Pi)Pi + (Pi,P3)^^3| - <y{^)E. 



(58) 



According to the response equations (pS)) . see also Eq. ((55t . the Lagrange parameters Pi,P2 are determined via the 
generalized linear Boltzmann equations, taken in Born approximation. 



Pi; Pi) ~iuj{Pi,Pi) 

uj-\-ie 

P3;Pi) -iio{P3,P,) 

LJ-\-ie 



Pi 
Pi 



P3;P3) -iw(Pi,P3) 
P3;P3) -iuj{P3,Ps) 

UJ-\-i€ 



P3 = (Pi,Pi)-/3P 



F3 = (P3,Pl)-/3P. 



(59) 



As shown in Appendix |D1 we have (A, A) = Nm/^, (A, A) = (A, A) = |iVm//3, (A, A) 



||7Vto//3 . Using 



Cramers rule, the response parameters Pi, P2 are expressed in terms of the electrical field E and correlation functions. 
For the dynamical conductivity, Eq. (1581) . we find after algebraic manipulations the expression Eq. (j57l) with 



zw7V^-f (Pi;Pi 



Pi; Pi 



Lj-\-ie 



P3;Pi 



P3;Pi 



(60) 



Evaluation of the correlation functions occurring in the renormalization factor r{io) in Born approximation is given 
in Appendix [El 



Results for the renormalization factor at solar core conditions and lower densities are shown in Fig. [TJ At 
solar core conditions (P 573 eV = 42.13 Ryd, n — 1.51 x 10^^ cm" -3 ^ 2.22 Gg^), we have a weakly inter- 
acting (plasma parameter F = e^/(47reo^BP) (47rn/3)^/^ = 0.1) and nearly degenerate (degeneration parameter 
Q = 2mkBT /h^ {3n'^n)~^^^ = 1.3) plasma. At the lower densities, the plasma becomes more classical. At high 
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FIG. 1: Frequency dependence of (a) the real part and (b) the imaginary part of the renormalization factor Eq. I|60p . Hydrogen 
plasmas at temperature T = 42.13 Ryd = 573 eV (solar core) and three different electron densities n are considered. 



frequencies (i.e. large compared with the inverse relaxation time), r{uj) approaches 1, and higher moments of the 
momentum distribution that describe the deformation from a shifted Fermi distribution are not relevant. In the 
static case, the real part Re r(0) shows the effect of e — e collisions according to the Spitzer result [1, H^. Since the 
Coulomb logarithm ([TTl) depends on the density, in addition to the correct prefactor also the density dependence of 
the Coulomb logarithm occurring in the different moments is seen. Only in the very low-density limit the different 
Coulomb logarithms cancel. 

So far we evaluated the equilibrium correlation functions occurring in the generalized linear Boltzmann equation 
Eq. pip with the help of perturbation theory. Thus we solved a kinetic equation using a variational approach or a 
reduced set of relevant observables. Note that one can go beyond the kinetic equation that treats the single-particle 
distribution function by considering fluctuations in the two-particle states as additional relevant observables in the 
generahzed LRT ^Eo]. 



D. Limiting cases 

1. Zero- frequency limit: Static conductivity 



We rewrite the dynamical collision frequency ()49|) in a symmetric form by transforming x 
expression and using the Dirac identity. 



in half of the 



Z43(w) 



.9clcgcn 



dy 



y 



dx ( V 
— < i:5{x — w) + 7r(5(.T + w) — i 

X \ X — 1 



V 



In 



1 + e^^^/y^y^^^^^^ 



The pricipal values compensate in the static case w = Q. After expanding for small x, e '^^/^ 
the integral over x can be performed with the result 



lim v-ci{lj) 



dy 



(61) 

-y'+/^M[i+2x], 



(62) 



Note that only e — i collisions contribute to the one-moment Born approximation. 

First we discuss the Lorentz model. It is solved for the static case in KT using an energy dependent relaxation time. 
The dc conductivity in Born approximation for the one-moment approach (j53p . crD(O) = eoWpj/i^nlO), is not identical 
with (Tdc obtained from Eq. (|T8)) with the Coulomb logarithm ([T7|) . because l/fD(0) ^ TLorontz- This difference stems 
from the fact that in the one moment approach with the variational parameter Fi the p dependence is specified 
as gi{p) = hpE- The p dependence necessary for the Lorentz model to reproduce the result for the relaxation 
time approach is given by gi{p), see Eq. (jS]), and is only roughly approximated by gi{p) within the interval of 
relevance. However, If we add further moments g^ip), not necessarily including g^ip), the approximation of the exact 
p dependence is improving. This has already been extensively investigated, see Refs. |37l. |41|). The dc conductivity 
within LRT follows from Eq. ([57| as (t(0) = €oU!^i/[r{0)iy'-^^\0)] with the static renormalization factor r(0). The 



12 



collision frequency i/*^-^i^(0) improves the Born approximation i^(0) if further effects like dynamical screening and 
strong collisions are included. 

The equivalence of the KT and LRT for the Lorentz plasma in the static case w = can be shown rigorously by 
inspection of the kinetic equation. Taking the linearized Boltzmann equation (|¥T|) with the collision term (1551) , ([5^ 
in the static limit, 

- ^ - fp)p- ^ - - E ^' (^)^P' (63) 

p' 

= IT E HE^t, - Bp) [Sf{p + f) - Sf(p)] (64) 

where the expression psp is used to insert the change of the single-particle distribution function Sf{p) after expanding 
fp+q — fp ~ {d / dj3Ep)fp = —j3{Ep+q — Ep)fp{l — fp). This equation coincides with the equation of motion for the 
single-particle distribution function (jlip . that is obtained in the static case from KT and is solved using the relaxation 
time ansatz. 

Considering the electron-ion plasma, it should be pointed out that the relaxation time approximation is not appli- 
cable if electron-electron collisions are relevant. In contrast, (t(0) obtained from LRT contains also the contribution of 
electron-electron collisions as given by Eq. (j40l) in the static limit. For this, the static renormalization factor rjO) can 
be evaluated from Eq. (|60p . In particular, it gives the correct Spitzer result if strong collisions are included [T9l[37ll4l|. 
see also Sec. IIVCI 



2. High-frequency limit: inverse bremsstrahlung absorption 

The dielectric function e{uj) — [^^(a;) + ic/(2a;)a(a;)]^/^ determines the refraction index nr{uj) as well as the 
absorption coefficient a{uj). We consider the long- wavelength limit where the transversal and longitudinal dielectric 
function coincide. The dielectric function or the optical conductivity cr(w) can be used to calculate the inverse 
bremsstrahlung absorption. In the high-frequency limit, where nr{uj) « 1 and w ^ i/, we have 

a{uj) = ^—^iTCiduj) « —^Jiev{uj) (65) 

so that the inverse bremsstrahlung absorption coefficient is directly related to the dynamical collision frequency 
obtained above from the solution of the Boltzmann equation. 

Bremsstrahlung radiation is described by the Bethe-Heitler expression resulting from QED in second order of 
interaction [4^|43j. In the non-relativistic limit and for soft photons, the absorption coefficient for a hydrogen plasma 
{Z, ^ 1) is given by gllii] 

where Ko{x) = dtexp[—x cosh(i)] = dy exp[— — x^/ {4:y^)]/y is the modified Bessel function of zeroth order. 

Generalized LRT gives the same result. We use the collision frequency Eq. ([5T|) in the nondegenerate case. At 
finite frequencies w, the integral with n = is no longer divergent at y = 0. Therefore, the screening of the Coulomb 
potential can be neglected {h = 0). We find P. 

with the free-free Gaunt factor in Born approximation 

= ^e^'^/'Ko . (68) 

The well-known Kramers formula for the inverse bremsstrahlung absorption [46!l results with the Gaunt factor 

^Kramers/, ,\ i 

9ff [^) - J-- 
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The one-moment Born approximation can be improved taking into account dynamical screening, strong collisions, 
and higher moments of the distribution function, as discussed earlier. However, in the high-frequency limit, the 
dynamical screening is not of relevance. The frequency dependence of the renormalization factor has been discussed 
in fJ], see also Fig. [TJ and converges to 1 in the high-frequency limit. Strong collisions have been considered and lead 
to the famous Sommerfeld result for the Gaunt factor [43,11^. For dense plasmas, the account of ion correlation S{q) 
[see Eq. ^] has a major effect and can directly included in the Born approximation j49| . 

The standard treatment of the kinetic equation using a relaxation time ansatz, see Subsec. Ill Bl fails to describe 
inverse bremsstrahlung absorption. The frequently used expression (|2ip for the dynamical conductivity, or the cor- 
responding expression for the dielectric function, are restricted to the low-frequency region since a static, but energy 
dependent relaxation time cannot be applied to the high-frequency region. Different approaches using Fermi's golden 
rule have been used [Tsj to derive expressions for the emission of radiation. A common treatment unifying both 
limiting cases, w — >■ and w — ^ c», is missing in KT within the relaxation time approximation. 

In contrast, our approach within LRT covers the entire frequency regime consistently. Note that it can also be 
applied to the degenerate case and to the relativistic regime, see [HO]- An important feature of the LRT is the 
possibility to include medium effects in dense plasmas such as the Landau-Pomeranchuk-Migdal effect [si'] . 



E. Dimensionless dynamical conductivity 

In the following we use Rydberg units where ti = 1, as = I, m = 1/2, e^/(47reo) — "2, ks = 1- The temperature 
T is then given in Ryd =13.6 eV and the electron density n in a^^. We introduce dimensionless quantities uj* = 
uj/ujp\ = wT/y/Trn and 

e2«3/2 1/2 

.*H = -|^.H. (69) 

In Fig. HJa), the ratio of the kinetic theory to the linear response theory is shown for the real part of the dynamical 
conductivity at various parameter values. The one-moment approximation is used, corresponding to the force-force 
correlation function. In Fig. El^b), the renormalization factor is included. In the low- frequency limit, deviations are 
shown that are due to the inclusion of e-e contributions. We give the limits of the expressions Eq. (IE5[) and Eq. (IE4|) . 
given in Appendix |E1 in the static case 

25/2 I 

CTkt(w = 0) 



7r3/2 Art 



<^LRT,l('- = 0) = 25/2^1/2 ■ ^^^^ 

In both approaches, the Coulomb logarithm behaves like lim„_j.o A ^ — ilnn in the low-density limit. At finite 
densities, different expressions are observed. The prefactor of the inverse Coulomb logarithm takes the value 1.015 for 
the Lorentz model that corresponds to the KT in relaxation time approximation. The Spitzer value 0.591 is approached 
in the LRT considering the Born approximation (0.2992 in the one-moment case, 0.5781 in the two-moment case). 
This quick convergence is known from the literature, see [sj]- The inclusion of the third moment of the momentum 
distribution takes electron-electron interaction as well as transport of heat into account. 
In the high frequency limit, we find from Eq. (jE5[) and Eq. (|E4[) the asymptotic expansions 



„ * / N l6V2n ^ 1 

RecrKT('^ ^ OO) = 3^y3 ^KT^ 

y2ni/4 1 

Rea£R,T,i(^^oo) = 3^5/4^3/2 ^7/2 ■ (^1) 

The ratio between KT and LRT behaves as uj^^'^. Thus, in the high frequency limit, the ratio diverges, see Fig. [51 In 
conclusion, above the plasma frequency the kinetic approach becomes essentially wrong. 



V. CONCLUSION 



Considering the interaction of radiation with matter, often a dielectric function or dynamical conductivity is used 
that is derived from kinetic theory using an energy dependent relaxation time, see Eq. (j2ip and Eq. (|E5|) . However, 
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FIG. 2: Ratio of the real part of the dynamical conductivity, calculated within relaxation time ansatz Eq. (|E5|) in comparison to 
generalized linear Boltzmann equation (LRT) in (a) one-moment approximation Eq. (|E4p and (b) two-moment approximation 
lEll Hydrogen plasmas at temperature T — 42.13 Ryd = 573 eV (solar core) and three different electron densities n are 
considered. 



this expression is valid only for elastic collisions of electrons so that electron-electron collisions cannot be included. 
Furthermore, the frequency dependence is not correctly described. In particular, the high-frequency behavior has a 
wrong dependence on uj and fails to describe inverse bremsstrahlung. We developed an alternative approach that is 
free from these shortcomings. 

We have derived a generalized linear Boltzmann equation Eq. (1311) that is valid for any frequencies and at arbitrary 
degeneracy. Besides electron-ion interaction, also electron-electron interaction is included. The drift term and the 
collision term are expressed in terms of equilibrium correlation functions that are, in general, complex quantities. In 
order to apply this approach consistently, one has to deal with two problems as follows. 

Firstly, the correlation functions can be evaluated numerically or, using quantum statistical methods, in perturbation 
theory. As simplest approximation, we considered the Born approximation, see Eq. (|53| with Eq. p9)) . Eq. ((5T|) . and 
Eq. (|Eip . This leads to analytic expressions that are tractable to be used for simple evaluations. 

Secondly, solving the generalized linear Boltzmann equation, a variational principle has been applied that optimizes 
the single-particle distribution function within a subspace of trial functions. In particular, we considered a finite 
number of moments of the distribution function. The single-moment treatment gives a result for the dynamical 
conductivity that is improved if higher moments of the distribution functions are taken into account. The contribution 
of higher moments is represented by the renormalization factor r{uj) that is, in general, a complex quantity. The 
high-frequency limit is not modified by the inclusion of higher moments and reproduces the well-known results for 
bremsstrahlung. The static limit converges to the Spitzer result for the conductivity with the inclusion of higher 
moments that desribe also the contribution of electron-electron interaction. 

We compared both approaches for different plasma properties. In the case of the Lorentz plasma that takes into 
account only elastic scattering of electrons by the ions, the correct static conductivity is obtained in KT using an 
energy dependent relaxation time. To get this result in LRT, the variational solution with only the lowest moment 
Pi is not sufficient, and higher moments should be considered. In particular, the inclusion of the fourth moment 
P4 alone gives the exact result for the static conductivity, the solution of KT with an energy dependent relaxation 
time becomes increasingly inappropriate with higher frequencies. In contrast, the expressions obtained from LRT are 
applicable at any frequency. 

Considering the more realistic case of the electron-ion plasma, the relaxation time ansatz to solve the kinetic 
equation breaks down. The inclusion of electron-electron collisions where the single-particle energy is not conserved 
represents no problem in LRT. The exact results for the transport coefficients in the low-density limit given by the 
Spitzer formula are reproduced by the LRT, in contrast to KT. The correct treatment of inverse bremsstrahlung shows 
that LRT is valid in the entire frequency domain, in contrast to KT using the energy dependent relaxation time that 
cannot reproduce the correct frequency dependence of the optical conductivity. 

Starting from a general LRT, a linearized Boltzmann kinetic equation has been obtained, and the relation to the 
results of th the relaxation time approach in the KT discussed. We restricted ourselves to a two-moment Born 
approximation. Possible improvements as pointed out throughout the paper are summarized here again as an outlook 
to further considerations and calculations. 
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• Taking the single-particle occupation number Up as relevant observables -B„, the deviations from equilibrium 
(rip)* — fo{p) describe the non-equilibrium state. The set of relevant observables can be extended by including 
initial state correlations jn particular the formation of bound states. This is straight forward in a general version 
of the LRT, see e.g. [illjls^] Sophisticated approaches have been worked out to show conservation of total energy 
and the systematic inclusion of correlations and bound state formation, using non-equilibrium Green's function 
theory [53, 113] or within generalized linear response theory [l^, |5^ . This is of relevance to investigate partially 
ionized plasmas, but also allows for the treatment of quasiparticle formation and the Debye-Onsager relaxation 
effect. 

• In linear response theory, the drift term and the collision term are expressed in terms of equilibrium correlation 
functions. They can be evaluated numerically or within perturbation theory, if we expand with respect to the 
interaction. The Born approximation is improved if higher orders with respect to the interaction are taken into 
account. The technique of thermodynamic Greens functions has been used for the evaluation of equilibrium 
correlation functions 0, The binary collision approximation is obtained if ladder diagrams are summed up. 
Dynamical screening results from the summation of ring diagrams. Perturbation expansions are more efficient 
if correlations are already included in the set of relevant observables so that they dan't have to be generated by 
a dynamical treatment, i.e., by considering higher order perturbation expansions. As example, we refer to the 
formation of bound states discussed above. Instead of finding their influence using higher orders of perturbation 
theory, we can treat them as new degrees of freedom introducing the corresponding relevant observables, e.g. 
their disribution function or a finite number of moments. Then, memory effects become less important, and the 
Markov approximation can be used, e.g., introducing stochastic forces [17^. 

Equilibrium correlation functions that determine the transport coefficients can be calculated for arbitrary frequen- 
cies, degeneracy, electron-electron collisions, and including collective excitations. The frequency dependence and 
further aspects are disregarded if a relaxation time is introduced. The relaxation time approach is exact only in 
the case of elastic scattering, for instance of electrons by ions in the adiabatic limit. Electron-electron scattering as 
well as finite frequencies of the electric field cannot be treated by the relaxation time ansatz. Thus, the generalized 
linear Boltzmann equation obtained from linear response theory reproduces some well-known benchmarks such as the 
Spitzer result for the static conductivity of the fully ionized plasma or the Kramers formula for the bremsstrahlung. 



Appendix A: Derivation of the response equations 



The hermitean observables Bn are assumed to conserve the total particle number so that the entropy operator 
H — fiN is replaced by the system's Hamiltonian H in the A dependence of the relevant statistical operator ([25]). Note 
that the averages are calculated with the equilibrium statistical operator that is known to us, and quantum statistical 
methods can be applied such as Green function techniques or numerical simulations to evaluate it. Thus, in linear 
response theory the Lagrange multipliers Fn{t) can be eliminated using equilibrium correlation fuctions. 

The relevant statistical operator serves as initial condition to determine the non-equilibrium statistical operator 
p{t). Further correlations are build up by the dynamical evolution [17] with the total Hamiltonian H^^^ = H + ^oxf 



m 



lim e 

e-i-O 



(Al) 



with the time evolution operator U{t,t') given by ih{d/dt)U{t,t') ~ Hl^JJ{t,t') and U{t,t) — 1. The external 
perturbation to the system's Hamiltonian H shall have the general form /f*xt — Decomposition of the 

time dependence of the field into Fourier components hj{t) = /ij(ti;)e~"^*/2 + c.c. = Re{/ij(a;)e~*'^*} is particularly 
convenient in linear response since the reaction of the system is the superposition of the reaction to different spectral 
components of the external perturbation. Subsequently, the time dependence of the response to each component will 
have the same frequency in the stationary case, i.e. Fn{t) = Re{i^„(a;)e~*"*}. In the following, we consider a fixed 
value w for the frequency of the external perturbation. 

We now perform a partial integration of the statistical operator (jAip and linearize with respect to the external 
fields hj and the response parameters Fn, 



Pirrel(t) = /5(t) - Pidit) 



lim 



f^i'e-«(*-*')e-»-^(*~*')/'i 
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H,6B„ 



Fr,{t')+5B^—F^{t') 



^l3X(H-^iN) 



(A2) 
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According to Eq. (|24|) we have Ti {Bn pi„ci{t)} = 0, for details see Finally, applying the Kubo identity 

/? / d\e'^'^"[C,H]e^^"po^ f d\^C{^in(3X)po ^ [C Jo] , (A3) 
JO Jo ""^ 

with C = i^cxt; fi^'i expression that relates the response parameters _F„ to the external fields hj, 

m ' j 

where the Laplace transform of the correlation functions (j29p has been introduced. After partial integration, 
—izlA-.Bj = {A,B) — IA]B\ = {A,B) + IA]B\ , we arrive at the response equations (|28l) with the exter- 



nal perturbation ff*^^ = —eR ■ E{t), R — r^, and R = P/m. 



Appendix B: Evaluation of the collision term 



We evaluate the Onsager coefficient Cpp'{uj) = {"hp' ; rip) ^^-^ which occurs in the collision term Eq. (|38p of the 
linearized equation of motion for the single-particle distribution function Eq. pip . Inserting the time derivative of 
the occupation number (using V*{~q) — V{q)), 



(Bl) 



•hp = ^[H,hp] ^ ^^Vci{q) 
q 

+ ^ ^°'=(^) [al+qO-l'-qap'O-p - alal,dp'-gap+q 

p'q 

into Eqs. (j27p . (j29p . we evaluate the correlation functions for the electron- ion contribution in Born approximation 

-2 VoiiqWM) / dt e^("+-'* / dA (B2) 
n Jo Jo 

I [Xr {poap+g«P«p'+g' v} - Tr |po 

r i^pQalap+gal,^^,dp, | - Tr jpoapOp+gap' } 



} 



QjiiEp-Ep+,){t-thp\)^^ . 



The A intergral can be executed. The application of the Wick theorem to the quantum statistical averages Tr {po- ■ ■} 
leads to 6 functions, in particular q = —q'. Contributions with q = cancel. We assume isotropic interaction 
V{q) = V{—q) and obtain 



1 ^ ^ef^iEp+.-Ep) _i 
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fi{Ep+q - Ep) 



i{uj + if) + i{Ep+q — Ep)!^ 
-1 



1? \ 1^'='^'^^' mv-Ep^q) ^"^^ ^'^''^W^) + ^{Ep - Ep+q)/h 



\Spf ^p-^q Sp/ ^p] 
\Sp' ^p Sp' ,p~\~q\ 



fp fp+q 

P{Ep+q - Ep) 



+ ie+ {Ep+q- Ep)/h u + ie- {Ep+q- Ep)/h 



i^p', p+q ^p' ,p] 



(B3) 



using (e^'-^p' — l)/p'(l — fp) = fp — fp'- Subsequently, the Onsager coefficient can be given as Eq. ([39)) . 
With this result, the collision term Eq. ([55]) for the Lorentz plasma reads 



P ^2 I 01 Wl ^(^^^^ -Ep)\uj + i€+ (Ep - Ep_q)/h ' co + ie- {Ep - Ep_q)/h 



{Fp+q - Fp) (B4) 
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which is now a frequency dependent and complex quantitiy. We can ehminate the Lagrange multipher Fp according 
Eq. psp in order to express the cohision integral in terms of the single-particle distribution function. 
A similar calculation gives the electron-electron contribution in Born approximation 

-fp+qfpi-q{'^ - /pi)(l - fp) 



P{Ep+q + Ep^^q — Ep^ — Ep) UJ + ie — {Ep+q + Ep^^q — Ep^ — Ep)/fl'' 



g/3(-Ep-(--Epi --Epi_5--Ep+5) _ 2^ ^ 

fi{Ep + Ep^ — Ep^^q — Ep+q) UJ + ie — {Ep + Ep^ — Ep^^q — Epj^q)/fi ^ ''^ ^^"^ 



X 



\5p\p+q + 5p' ,pi—q ^p',Pi ^P',p] I (-B5) 



where 14c,ox(9;p,Pi) = Vcc{q) ~ (Jcr^.o-sKicdpi — p — (fl) is the exchange interaction with ai denoting the spin explicitely. 
The respective Onsager coefficient can be given as Eq. (|40p . It is easily seen from the final expressions (pQ)) , (j40|) that 
the real part of the Onsager coefficient £.ppi{uj) — £°p, (w) + £pp,{(jj) is non-negative, ReCpp'{uj) > 0. 

Appendix C: Proof of the variational solution 

To begin with, we show that the entropy production Eq. (pSj). 

SintlGp] = J2G;iCpp,iuj) + c;,piuj))Gp, = ^g; v)^^^^ Gp, (ci) 
pp' pp' 

as a functional of an arbitrary Gp is positive definite. Using the spectral density of the operator G — ^pGphp, we 
find 

^int [GJ = (G;G) ^ ^s(Lu + -{E,,-E„,,))\{n\G\m)\^>0. (C2) 

Now we consider the functional Eq. (jCll) for the function {Gp — Fp) and decompose 

Sint[{Gp — Fp)] = Sint[Gp] — [G*Cpp' {uj)Fpi + c.c. ] - ^[F;/:pp.(c^)Gp, + c.c] + Sint[Fp] (C3) 

pp' pp' 

Making use of the constraint Eq. the first contribution is expressed as 

^int [Gp] [g; + Gp] Dp, (C4) 

p 

the terms with iQp compensate. Since Fp solves the linear Boltzmann equation (|42p . the second contribution is 
transformed into 

J2 G*pCpp' {uj)Fp. + c.c. = ^pDp + Y i^pG*pFp + c.c. (C5) 
pp' p p 

For the transformation of the third term, we use the symmetry Cpp'{uj) = Cpip{uj) due to detailed balance which 
can be seen easily from the explicit expressions Eq. (|B3I) and Eq. (IBSp . Furthermore, the proof of the reciprocitivity 
condition Cpp'{uj) — Cpp,{—uj) can be shown generally using the eigenstates \n) of the system Hamiltonian, 

interchanging n and m. Finally, we find 

5]F;/:pp-(c.)Gp, = ^[G;,£;p,(c.)Fp]* = 5][G;£pp,(~u;)^pT = 5^[Dp -f zf7pF;]Gp . (C7) 

pp' pp' pp' p 
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We sum up all contributions in Eq. (103)) using the Eqs. ((C4)) . ((C5|) . (ICT)) . 

^i„t [{Gp-Fp)]^ 5i„t [fp] - 5i„t [Gp] > . (C8) 



This is a positive definite expression due to Eq. (jC2[) . Thus we find that the entropy production is maximal if the 
trial function Gp is the solution Fp of the Boltzmann equation. 

Appendix D: Evaluation of Eq. (|47|) 

We execute the p integration on the left hand side of Eq. (|47l) with — 

iV?iVffl-fl ^^J}o_ fj^ ( dfp \ , _ Anm (2m)3/2 f 3/2 

3 2.^P/p(i /pj - 3 (2^)3 y^fl^ cJ/Si^JP'^^- 3/3 (2^f j dEp^P "'^^ 



~ a^(2^ 



after integration by parts. This is also identical with (^Pi,Pij which is the Kubo scalar product (l?7l) of the first 
moment ([5]). 

In the collision term, that is the second term on the right hand side of Eq. (ITT)) , we insert the expression Eq. p9l) . 
The sum over p' is immediately executed and gives qE- The first contribution ( from 5 function) as well as to the third 
contribution (from first principal part) are considered together and can be transformed hy q ~q, then p ^ p + q, 
so that they coincide with the second and fourth contributions, respectively. We find after canceling some common 
factors 



From Eq. (HH) we find 



eh^NE - ~tuj-n^NFi - V |K=i(<z)|'g| V . ^^^i • (D2) 



mN ^ P[Ep+q - Ep) LJ + ie+ [Ep+q ~ Ep)/n 

We shift p — 7- p— q/2 so that Ep^g/2 — -E'p-q/2 = • q/fn and with spin factor 2, 

mA ^ ^^i^g {^T'j J -00 s uj + ie + Tiqs/m 

l-oo g/3(?iV2m)(r2+s2+sg+,2/4)_/3^ _ ^l3(h^ /2m){r'^ + s'^ -sq+q^ /4)-Pti 







where cylindrical coordinates with respect to the q direction have been introduced, s is the component of p in (f 
direction, r is the component orthogonal to this axis. The integral over r can be performed, 

I r 1 1 

ar 



2 _/q g/3(ri2/2)n)(r2 + s2+sg+g2/4)_^^ ^ -j^ g-/3(R2/2m)(r2+s2-sg+g2/4)+/3/i ^ ]^ 

^ 1 ^ _|_ g-/3(nV2'n)(s-g/2)^+/3Ai" 

■ In 



^7j2 ^p{h^/m)sq _ I 



I ^ g-/3(fi2/2m)(s+g/2)2+;3M 



(D5) 



Furthermore, we neglect the ion correlation so that S{q) — 1 for the structure factor. Note, that the Born approx- 
imation, Eq. (|D4I) , is divergent at zero frequency. As well known, this problem is solved if we go beyond the Born 
approximation and take into account higher order contributions due to dynamical screening and strong collisions. 
This was already discussed by Landau and Lifshitz [56j and has been shown to be consistent using Green function 
techniques, see [3ll3.[29j. In this way, the correct zero frequency limit of the collision frequency is obtained. As shown 
in Ref. [13], alternatively, the Coulomb potential in Eq. (jD4p can be replaced by a statically screened potential, the 



Debye potential dH]), so that |Voi(g)P - NV^. With s = Jj^^.q= J^V. expression (03) follows 
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Appendix E: Renormalization factor and dynamical conductivity 

We use Rydberg units as introduced at the beginning of Section HVE] and in Eq. §9^1 . In LRT, the conductivity [57] 
within one-moment Born approximation in the non-degenerate hmit [5T] gives (w — u* ^T:n/T) 



2n 



VjJ 



dy 



(y2 + 27rn/T2)2 



dx- 



-4:Xy 



xy{w — xy + ie) 



(El) 



The renormahzation factor r(w) is taken with the first and third moment of the distribution function (i.e. particle 
current and energy current). According to Eq. (j60p . generahzed force-force correlation functions have to be calculated 
after decomposition: {Pu Pm)uj+it — {FT'^ Prnli^+i^ + {Pi'^'j Pm)'-i+i<^ ■ Considering non-degenerate limit of the Born 
approximation again, we have from the electron-ion interaction 



/ pci. pen 



4 Nn 



dy 



^V^VtJo " (y2 + 27rn/T2)2 



dx 



-Axy 



xy{w — xy + ie) 



{x,y} 



Im 



(E2) 



where {x, y}";^ = 1, {x, yjg^ = l + 'ix^ +y^, and {x, yY^^ = 2 + IQx^ + 9a;'* + 2y^ + Qx^y^ -f y"^ 
For the electron-electron interaction we find 



/ pcc . p. 



Nn 



rn /uj+ie 



3V2tt VT 



dy 



y 



{y^ + Aim/T^y 



dx- 



1 



-Axy 



xy{w — xy + ie) 



-ix-yf 



{x,y} 



Im 



(E3) 



where due to momentum conservation {Pi'^ — 0) we have {0;,?/}°'^ — {x,yY^-^ — and {x.yjgg = 1 + (19/4)a;2 



For the evaluation we use 



-|- i7rS{xy — w). The 6 function allows to perform the integral over x 



xy—w—ie xy — w 

to obtain the real part of the correlation functions {Pf, Pm)u)+if For the imaginary part, we also can perform the x 

integral after partial fraction decomposition and using V J^^dx^ = Tie-'' er?i{a). 
In particular, we have for the single moment approximation where r{w) — 1 



'LRT,1 



2n 



iy' 



27m \2 



2i e 



-(y- 



erfi(y 



w 

- - e 

y 



erfi(j/) 



(E4) 

-1 



For direct comparison, we give explicitely the dynamical conductivity from KT (|2ip with the energy dependent 
relaxation time for the Lorentz plasma (jl7p 



8 



^KT 



2n 1 



dx 



30F V 7rT3 T 7o iw* - [ln(l + &) - 6/(1 + ^)] 1^^ 



(E5) 



with h = x'^T/{2TTn). 
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